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Abstract 

Entropy production is one of the most important characteristics of non-equilibrium steady states. 
We study here the steady-state entropy production, both at short times as well as in the long-time 
limit, of two important classes of non-equilibrium systems: transport systems and reaction-diffusion 
systems. The usefulness of the mean entropy production rate and of the large deviation function 
of the entropy production for characterizing non-equilibrium steady states of interacting many- 
body systems is discussed. We show that the large deviation function displays a kink-like feature 
at zero entropy production that is similar to that observed for a single particle driven along a 
periodic potential. This kink is a direct consequence of the detailed fluctuation theorem fulfilled 
by the probability distribution of the entropy production and is therefore a generic feature of the 
corresponding large deviation function. 

PACS numbers: 05.40.-a,05.70.Ln,05.20.-y 
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I. INTRODUCTION 



The study of large deviation functions has been of increasing importance for the un- 
derstanding of many-body systems. On the one hand large deviation functions form the 
basis of a modern approach to equilibrium statistical mechanics [l-3|, on the other hand 
they are increasingly recognized of being of fundamental interest for the characterization 
of non-equilibrium systems where they are sometimes considered to play a role similar to 
that played by the free energy at equilibrium. Important examples are given by the large 
deviations of the steady-state currents 0, 0] in systems that are far from equilibrium. 

Recently, Mehl et al. [oj] extended the study of large deviation functions far from equilib- 
rium to the steady-state entropy production. Studying a one-dimensional system composed 
of a single particle driven along a periodic potential, they reformulated the problem as a 
time-independent eigenvalue problem j3| and showed that for this system the large devia- 
tion function of the entropy production exhibits a kink at zero entropy production. Similar 



arge deviation 
4l3H. These are 



kinks in the large deviation function of the entropy production or in related 
functions of other quantities can also be found in a range of other systems 
intriguing results that raise the important question whether the presence of this kink is a 
universal feature of systems with non-equilibrium steady states. 

In this paper we study the entropy production in the steady states of different non- 
equilibHum interacting „.any-body =yste,ns^On the one hand we =tudy the Patt.ally A^ym- 
metric Simple Exclusion Process (PASEP) [14], a transport process in an open system where 
particles can enter or leave only at the system boundaries. On the other hand we investigate 
react ion- diffusion systems on a ring where the number of particles can change everywhere 
in the system due to reactions between the particles. Our main emphasis thereby is on the 
entropy production in the long-time limit where we study the large deviation function in 
the same way as done by Mehl et al. in their study of the single particle system. For all 
studied many-body systems we find a kink-like feature at zero entropy production, similar 
to what has been observed for the single particle system. We show that this kink is a generic 
feature in non-equilibrium systems obeying a detailed fluctuation theorem for the entropy 
production. 

It is worth noting that large deviation functions of the current have been studied previ- 
ously in some related systems. Studies of the Totally Asymmetric Simple Exclusion Process 
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(TASEP) [15|, ll6|] and of the Symmetric Simple Exclusion Process (SSEP) revealed a 

non-analytical behavior of the current large deviation function at the value of the parame- 
ter that is canonically conjugated to the particle current. Using dynamical renormalization 
group techniques, it was later shown that the value of the corresponding exponent follows 



from the noise renormalization [18]. Very recently, Bodineau and Lagouge 19| studied the 



current large deviations in a driven dissipative lattice gas model with creation and annihi- 



lation processes, whereas Simon 20| investigated the large deviation function of the current 
for the Weakly Asymmetric Simple Exclusion Process on a ring. 

Our paper is organized in the following way. In the next Section we introduce our models, 
before discussing in Section III our approaches for computing the entropy production at short 
times as well as the mean entropy production rate and the large deviation function for the 
entropy production in the long-time limit. Our results are presented and discussed in Section 
IV. Finally, Section V gives our conclusions. 

II. MODELS 

In order to elucidate the entropy production in interacting many-body systems, we discuss 
in the following reaction-diffusion systems on a ring as well as open transport systems where 
particles move through a system that they can enter or leave only at its boundaries. In the 
past, due to the combination of their conceptual simplicity and highly non-trivial results. 



diffusion- limited reaction systems [2l|,|22| and simple exclusion processes |23l-l25| have greatly 
contributed to our understanding of processes far from equilibrium. For the same reasons, 
these models are also the natural choices for our study of the entropy production in many- 
body systems. 

As react ion- diffusion systems we consider simple cases where particles A diffuse on a one- 
dimensional ring, under the condition that every lattice site can only be simply occupied. 
This diffusion process is mimicked by the hopping of a particle to an empty nearest neighbor 
site with rate D. In addition, we also allow for particle creation and annihilation. In the 
creation process a new particle can be created at an empty site with rate h, whereas in the 
annihilation process n particles on n connected sites are destroyed with rate A, yielding n 
empty sites. In the following we characterize our models by the number of particles involved 
in the annihilation process and call M„ (with n = 2,3) the model where n particles are 
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destroyed. As a variant, we also study the situation (we call the resulting model M2/) where 
in a two particles reaction only one particle is destro yed . 
These models are the same as those studied in 
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28l ] in order to better understand 



steady-state and transient properties of reaction-diffusion systems. All these systems have 
non-equilibrium steady states. In addition, some or all reactions do not allow for a direct 
back-reaction, and microscopic reversibility is broken. However, as discussed in the next 
Section, the computation of the entropy production requires microscopic reversibility, i.e. 
for every reaction the direct back-reaction must be possible. For that reason we are adding 
to the reaction schemes the direct back-reactions of the reactions just described. Thus, any 
particle can be destroyed with rate Eh h, with < e/j < 1, irrespective of whether neighboring 
sites are occupied or not, and a reversing of the annihilation process yields the creation of 
particles with rate A, where < < 1- It is important to note that even with these 
expanded reaction schemes our systems are still characterized by non-equilibrium steady 
states, as long as eh,£x < 1. In the following we set Eh = £x = £ and focus on the case 
£ <^ 1. For the convenience of the reader we summarize the different reaction schemes in 
Table E 



M2' 




A 

A + A ^ + A 

h 

0^ A 

£hh 


A 

nA ^ nO 

e\X 
h 

Qt± A 

Ehh 



TABLE I: The different reaction schemes discussed in this work. The back-reactions are taking 
place with rates Shh and ex\, with Q < Sh < 1 and Q < e\ < 1. 



As transport process we consider the Partially Asymmetric Simple Exclusion Process 
(PASEP) where every configuration is reversible. In the Total Asymmetric Simple Exclusion 
Process (TASEP) in an open one-dimensional system, particles that are fed into the system 
at, say, the left end with rate a can leave the system at the right end with rate /3. Inside the 
system particles can only jump to a right neighboring site, provided that site is not occupied. 
Obviously, microscopic reversibility is always broken for this model. In the PASEP, however, 
particles can always jump in both directions, provided that the chosen site is empty, with 
the probability for a jump to the right being p, whereas the probability for a jump to the 
left is 1 — p. In addition, particles can exit the system at the left with probability SaCx (with 
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< ^Q, < 1) and enter at the right with probabihty (with < Ejs < 1). It follows that 
microscopic reversibility is always fulfilled. In order to keep the number of parameters as 
small as possible, we set Sa = £/3 = £, with e 1. 



III. METHODS 



Whereas our main focus in the following will be on the large deviation function of the 
steady-state entropy production in the long time limit, we will also briefiy discuss the entropy 
production at short times. 

Focusing on our lattice models, let us consider a path in configuration space Cq — > 
Ci — )■ ■ ■ • — )■ Cm-1 — > Cm that starts at some configuration Cq and ends at some 
configuration Cm after M diffusion or reaction steps. Every configuration Ci is uniquely 
characterized by the occupation numbers (being or 1) of all lattice sites. We will denote 
by Ps{Ci) the probability to find the configuration Ci in the steady state. With every step 
i leading from configuration Cj_i to configuration Ci we associate a time increment given 
by 

(1) 



where u{Ci^i — > Cj) is the rate with which we go from configuration Cj_i to any other 
accessible configuration Cj. The sum in the denominator is thereby a sum over all configu- 
rations Cj (including the configuration Ci in which the system will be at step i) that can be 
reached from the configuration Cj_i through diffusion or reaction. With that we can assign 
a total time 

M 

1=1 

to our trajectory in configuration space. 



Along the same trajectory the total entropy production is given by 29| 

= lu|f°\ + .,„ (3) 

where is the entropy produced in the particle bath connected to our system. The bound- 
ary term, In-^^^^, which can be neglected in the long time limit, needs to be included 
when investigating fiuctuation relations of the entropy production at short times. 



It is worth pointing out that in Eq. (|3]) the transition rates from configuration Cj_i to 
Ci and from configuration Ci to Cj-i are showing up together. This is the reason why we 
have added the back-reactions to our reaction schemes. 

In Section IV. A we will discuss the mean entropy production in the short time regime 
as a function of time r. For very small r, i.e. values of r where only few configurations 
are visited, this can be done in an exact way, where, starting from all possible initial states, 
all possible trajectories are determined and the corresponding entropy changes are recorded 
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2l\ . This allows us to determine in a numerically exact way the probability distribution 
of the total entropy production, p{stot, after time r. For larger (but still small) values of 
r, we can compute the mean entropy production through standard Monte Carlo sampling. 
All the data discussed in the following are numerically exact data, with the exception of the 
data shown in Fig. [2] which have been obtained by Monte Carlo simulations. 

With the definition ([3]) of the total entropy production the probability distribution 
p{stot, t) obeys a detailed fiuctuation theorem for any time r: 

p(stot,i")/p(-Stot,r) = exp (-Stot) . (4) 

Initially [tI. I30I32I this relation was shown to be valid in the long-time limit when only Sm is 



considered. If one adds the boundary term, as done in Eq. ([3]), then the theorem (jlj) holds 
true for any finite time interval r 



In the following much emphasis will be put on the large deviation function of the entropy 
production (also called the rate function) 



X(a-) = lim 



T— >00 



--lnp(sm,r) 
r 



(5) 



which describes the asymptotic large fiuctuations of the entropy production (in that limit 
we can neglect the boundary term in ([3])): 

P(sm, r) ~ exp (-x(CT)r) (6) 

with the normalized entropy production rate 

(y = Sm/ (r (s„)) . (7) 

Here (s^) = (^f") is the mean entropy production rate. In order to compute the large 
deviation function, we proceed as in Ref. [7| and consider the generating function 

(exp [-fiSmir)]) = j dsm exp {-^iSm{r)\ p{sm, t) (8) 
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of Sm(T). The time evolution of the generating function is given by the hnear operator 



—/I In 



^ Ci)e -r{Ci) 

j 



(9) 



where r{Ci) = J^^i^i — ^ ^j)- The spectrum of is strictly positive for yU 7^ 0, whereas 

j 

for /i = the lowest eigenvalue is zero and characterizes the steady state of the system. In 
the long time limit the lowest eigenvalue z/(/i) can equivalently be determined as 



lim 

r— ^oo 



--In (exp [-/is„(r)]) 



T 



(10) 



For our finite systems we determine the eigenvalue i^(/i) through the numerically obtained 
spectrum of the operator which is expressed by a finite matrix. This is done using 
standard algorithms. 

This lowest eigenvalue allows us to obtain expressions for various quantities of interest. 
Thus the mean entropy production rate is given by 



whereas the large deviation function is given by the Legendre transform 



X(cr) = max{z/(/i) - (sm)cr/i}. 



(11) 



(12) 



The fluctuation theorem 0, Isol i 
deviation function, 



32 1 is revealed through a symmetry relation for the large 



or, on the level of the eigenvalue, through the equivalent symmetry 

i/(/i) = 1/(1 - fi) . 



(13) 



(14) 



IV. ENTROPY PRODUCTION 

A. Entropy production at short times 

Before discussing the long-time behavior of the entropy production, which is the main 
focus of the paper, we briefly discuss some properties of the entropy production in the short 

7 



time limit. As an example, we discuss the reaction-diffusion models, but corresponding 
results are obtained for the transport processes. 

For all studied cases the probability distribution p{stot) for the total entropy production 
at short times is very irregular, as shown in Fig. [T^ for the model M3. Similar irregular 
shapes, dominated by prominent peaks, also show up when computing a related quantity, the 
so-called driving entropy production, in systems that are driven out of stationarity through a 



change of system parameters [26 
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34| . The details of the probability distributions depend 



on the microscopic details of the studied systems, especially on the values of the diffusion 
and reaction rates. Still, even though the probability distributions are very irregular and 
specific for the studied systems, we recover in all cases the celebrated fiuctuation theorem 
dl]), see Fig. [T]d, provided that the boundary terms are included in the definition of Stot- 

20 




FIG. 1: (a) Probability distribution for the total entropy production and (b) corresponding fluc- 
tuation theorem for model M3 on a lattice witli N = 8 sites. The values of the parameters are 
M = 6, D = 1, e = 0.01, A = 1, and /i = 0.5. 



With the probability distribution of the total entropy production at hand, we can straight- 
forwardly compute the mean value of entropy change (stot) during some time interval r. As 
shown for some cases in Fig. |2^, the total entropy production rapidly yields a linear increase 
with time. We also note, see Fig. |2]d, that the entropy change per site already becomes 
independent of the system size for very small systems. As a consequence of these two obser- 
vations we can restrict ourselves to small systems when discussing the long-time properties 
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T N 

FIG. 2: Mean entropy production in the steady state as a function of the total measurement time 
r (left) and system size N (right) for models (circles), M2 (squares) and M3 (diamonds). The 
parameters are D = 1, h = 0.5, A = 1 and e = 0.01, with (a) = 8 and (b) r = 5 for the right 
panel. Both the average entropy production per unit time, see inset in (a), as well as the average 
entropy production per site reach rapidly a constant value. 

in the following. 

B. The long time limit 

As discussed in Section III the long-time properties of the entropy production are encoded 
in the eigenvalue i/(/i) of the linear operator L^. After having determined i/(/x) we can use 
that quantity in order to discuss the mean entropy production rate, see Eq. ffTTj) . as well as 
the large deviation function for the entropy production, see Eq. 0121) . 

1. The mean entropy production rate 

In Fig. [3] we show examples of changes of the mean entropy production rate (sm) for 
various systems when varying one of the system parameters. Fig. [3^ illustrates some of our 
results for the react ion- diffusion systems where the creation rate h is changed and all the 
other system parameters are kept constant. Fig. [3)d is devoted to the transport process with 
a varying entrance rate a. 
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Interestingly, (sm) encodes important information on the physical properties of the differ- 
ent steady states. Focusing first on the reaction-diffusion systems, we see that for model M2 
the mean entropy production rate has a maximum around h = 1 and decreases towards zero 
for larger values of the creation rate. In stark contrast to this (sm) increases monotonously 
for the other two models M2 and M3. In order to understand this difference in behavior, 
we refer the reader to the reaction schemes listed in Table [B Noting that for very large 
values of h there is an increased probability for a free site to be surrounded by occupied 
neighboring sites, it follows that the creation process — )■ A is then effectively identical to 
the process A — )■ 2A. This is, however, exactly the reversed reaction to the annihilation 
process of model M2, such that this system approaches an equilibrium system for very large 
values of h. Consequently, the mean entropy production rate vanishes in that limit as no 
entropy production takes place in equilibrium. This is different for models M2 and M3, as 
here the reaction A — )■ 2A is not the direct back-reaction of the creation process, and the 
systems remain out of equilibrium for any value of h. 

Similar conclusions are also drawn when varying other system parameters. For example, 
when e — > 1, all reactions become symmetric, and the system approaches equilibrium. 
Concomitantly, the mean entropy production rate approaches zero in that limit. 

It is worth noting that the behavior of the mean entropy production rate is similar to 
that of a quantity derived from the stationary probability currents 35rl37| that has been 
proposed as a possible measure for the non-equilibrium character of a system. For example, 
both quantities show a non-monotonous change for model when increasing h, whereas 
their increase is monotonous for models M2 and M3 (see Fig. 2 in 26|). 

Non-trivial information is also contained in the mean entropy production rate of the 
transport process. As an example we show in Fig. [3)d (sm) as a function of the entrance rate 
a for a strongly asymmetric case with p = 0.95. The striking feature in this figure is the 
appearance of kinks in the curves obtained for smaller exit rates, as for example /3 = 0.25 
or /3 = 0.30. These kinks, which are located at a = /3, yield a sharp peak in the negative of 
the second derivative of (sm) with respect to a, see inset of Fig. [3b. We checked that these 
kinks show up for all /3 ^ 0.46 at the value a = (5. For larger values of /3 the kink vanishes 
and only a small peak remains in the second derivative (the peak still exists for /3 = 0.75 but 
is not visible on the scale of the inset in Fig. [3]d). Interestingly, in this regime the position 
of the peak is independent of the value of /3 and is given by a ~ 0.47. 
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FIG. 3: (Color online) Mean entropy production rate (sm) for (a) reaction-diffusion systems M2 
(circles), M2 (squares) and M3 (diamonds) as a function of the creation rate h (values of other 
parameters: N = 8, D = 0.5, e = 0.01, with A = 5 for model and A = 0.1 for models M2 and 
M3), and (b) the PASEP as a function of the entrance rate a (full black line: exit rate f3 = 0.25; 
gray line: /3 = 0.30; dashed black line: (3 = 0.5; green (dark gray) line: (3 = 0.75) with = 8, 
p = 0.95 and e = 0.01. Inset in (b): the negative of the second derivative of {sm) with respect to 
a reveals the existence of a kink in the mean entropy production rate for a = /3 ^ 0.46. For larger 
values of f3 only a small peak remains, the position of which is a ~ 0.47, independent of the value 
of /3. 



In order to understand this result we need to remind the reader of the phase diagram of 
the PASEP 



3^- For our choice of entrance and exit rates at the boundaries of the system 
the phase diagram of the PASEP is very similar to that of the TASEP and contains a low 
density phase, a high density phase and a maximum current phase. The phase boundary 
between the low and the high density phases, which is discontinuous and characterized by 
the coexistence of these two phases, is given by a = /3, a < Vc, where rc depends on the 
specific values of the different rates. In addition, one has that the boundary between the low 
density (high density) phase and the maximum current phase is given by a = and /3 > 
{a > Tc and (3 = rc). Insertion of the values of the rates used in Fig. [3b into the exact 
expressions for the phase boundaries 3^ yields the value Vc = 0.4545. Comparing this with 
the position of the kink in the mean entropy production rate for /3 ^ 0.46, we see that the 
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kink coincides with the discontinuous transition separating the low from the high density 
phase. In addition, for /3 > Tc the extremum of the second derivative of (sm) close to a ~ Tc 
indicates the transition separating the low density phase from the maximum current phase. 



2. Large deviation function for the entropy production 

In order to compute the large deviation function for the entropy production, we proceed 
as discussed in Section III and determine the lowest eigenvalue z/(/i) of the linear operator L^, 
see Eq. The large deviation function is then obtained through the Legendre transform 

Typical results obtained for our systems are shown in Fig. H] where the left column is 
devoted to the reaction-diffusion systems and the right one to the PASEP. As expected, the 
eigenvalue fulfills the relation (|T4l) that follows directly from the fluctuation theorem, 
see Fig. HIi and Fig. Hli. By construction, the large deviation function x((t) is zero at a = 1. 
In addition, plots of xi^^): see Fig. |Dd and Fig. H^, reveal for all studied situations a kink 
at zero entropy production cr = 0. This kink is best revealed when studying the derivatives 
of X with respect to a, as shown in Fig. and Fig. H^, yielding a step at a = in the first 
derivative and a maximum at cr = in the second derivative. 

As stressed in the Introduction, a kink has also been observed recently in a study of the 
steady-state entropy production for a single particle driven along a periodic potential 6|. 
The fact that a similar behavior is observed in our interacting many-body systems suggests 
that this feature is a universal one. 

Before discussing the origin of the kink further, let us briefly pause and discuss how the 
large deviation function and its derivative depend on the parameters of the system. As an 
example, we show in Fig. Elx and x' fo^' the react ion- diffusion model M3 for various values 
of the creation field h. In our discussion of the mean entropy production rate we pointed 
out that for this system (sm) increases when h increases, resulting in an enhanced non- 
equilibrium character when taking (sm) as a measure. Interestingly, the kink in the large 
deviation function also gets more pronounced when h increases, as witnessed by a larger 
amplitude of the step in the derivative x'- This is the generally observed behavior for all 
our systems: when the system parameters are varied in such a way that the mean entropy 
production rate increases (decreases), then the step height at cr = of the derivative of the 
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FIG. 4: (Color online) (a,d) The lowest eigenvalue vi^ii) of the linear operator L^, (b,e) the large 
deviation function for the energy production x^'^\ and (c,f) the derivative of the large deviation 
function for (a,b,c) the reaction-diffusion models shown in Fig. [3^, with N = 8, h = 2, D = 0.5, 
e = 0.01, A = 5 for model M2 and A = 0.1 for models M2 and M3, and (d,e,f) the transport 
processes discussed in Fig. [8)3, with N = 8, a = 0.6, p = 0.95 and e = 0.01. In all studied cases, 
the large deviation function exhibits a kink at zero entropy production, a = 

large deviation function also increases (decreases). 

In fact, the kink at a = in the large deviation function for the entropy production 
is a generic feature and immediately follows from the fluctuation theorem % |30|-|32|. To 
see this, we remind the reader that the fluctuation theorem shows up on the level of the 
large deviation function through the symmetry relation ( 1T3|) . Taking the derivative of that 
relation with respect to a yields the following relations between the derivatives at both sides 
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FIG. 5: (Color online) (a) The large deviation function of the entropy production and (b) its 
derivative for model M3 and different values of the creation rate h. The parameters used for this 
figure are N = 8, D = 0.5, A = 0.1, and e = 0.01. For cr = a kink shows up in the large deviation 
function which is best illustrated by the rapid change in its derivative with respect to a. 



of a = 0: 



da 

dx 
da 



(T=(70 



da 

dx 
da 



-O-Q 



- {Sm) 



2 — 
da 



(15) 
(16) 



CT=(TO 



valid for any value of the entropy production a^. In Fig. Elwe illustrate for model M3 the 
relation f[T51) for the large deviation function as well as the relations (fT^ and ffTBl) for its 
derivative. 

The relations (fT5|) and (fT6l) completely characterize the kink. In fact, one really should 
describe the observed feature as kink-like, as the derivative of x is still continuous at a = 
even though its value changes rather dramatically when the entropy production approaches 
zero. Eqs. f[T^ and f[TB|) make obvious the connection between the kink and the mean 
entropy production rate. It is the value of (sm) that characterizes this intriguing behavior 
of the large deviation function of the entropy production. A larger value of (sm) yields a 
more pronounced feature and a sharper and higher peak in the second derivative of the large 
deviation function. 

One might object that in some known cases the distribution of the entropy production 
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FIG. 6: (Color online) (a) Illustration of the symmetry relation (jl3p for model M3 with = 8, 
D = 0.5, e = 0.01, h = 2, and A = 0.1. The full line is ^(o')) the dashed line is x(— cr), whereas 
the green (gray) circles follow from x(~^) ~ <y{sm)- For the chosen set of parameters we have 
that (sm) = 11.1228. (b) Illustration of relations (jlSp and (jl6p of the derivative x' of the large 
deviation function. Full black line: x'('^); dashed black line: \'{(t) +x'{—(t); green (gray) line: 
x'('7) - x'(-f'"); black circles: 2\'{a) + {sm)- 

has a nearly Gaussian form and that for these kink-like behavior is not observed. 

An example can be found in Ref. js] (see their Figures 1 and 3), where the large deviation 
function for a particle driven along a periodic potential is nicely fitted to a parabola 

x(a) = ((U/4)(^-lf (17) 
for small forces. In a similar way, the largest eigenvalue seems to be well described by 

= (sm)/i(l - /Lt) • (18) 

These results seem to indicate that close to equilibrium systems with a nearly Gaussian 
entropy production can fulfill the fluctuation theorem without exhibiting a kink-like feature 
in the large deviation function. 

Our systems allow us to address this point in a systematic way. In fact, we always find 
for systems that are close to equilibrium (as characterized by a small value of the mean 
entropy production rate) that the entropy production approaches a nearly Gaussian form 
and that the kink seems to vanish. This is illustrated in Fig. [7] for the transport process with 
a = (3 = p = 0.5 and various values of e. In this case the system approaches equilibrium 
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FIG. 7: (Color online) Discussion of the transport process with a = (3 = 0.5 and p = 0.5. The 
different curves correspond to various values of e: 0.001 (dashed black line), 0.01 (dashed gray 
[magenta] line), 0.1 (full black line), and 0.5 (full gray [magenta] line), (a) Lowest eigenvalue 
z^(/i) of the linear operator L^, (b) large deviation function for the entropy function x{'^)^ (c) first 
derivative of the large deviation function, (d) second derivative of the large deviation function, 
and (e) second derivative of the large deviation function for e = 0.5. The dotted lines result from 
fitting the data for e = 0.5 to the quadratic forms (fT7|) and ([H]). Systematic deviations between 
the data and the fit are readily observed in }("{(t). 

when e increases, yielding an equilibrium system for e = 1. As we see from Fig. [7^ and 
[7)d, the kink-like feature in the large deviation function gets less and less pronounced when 
increasing e and seems to have vanished completely for e = 0.5. Indeed, in that case one 
can fit the largest eigenvalue and the large deviation function to the quadratic forms ( !T8|) 
and ( fT71) . thereby obtaining an excellent agreement. However, a closer look at the second 
derivative, see Fig. [711 and [7^, reveals that qualitatively the case e = 0.5 does not differ from 
the other cases: the second derivative exhibits a maximum at a = which is the signature 
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of the kink-like feature. Admittedly, this maximum is rather shallow for a system close to 
equilibrium, but it persists as long as the mean entropy production rate does not vanish. 
Obviously, it is not good enough to simply fit the data to the quadratic forms (IT7I) and 
f[T5]) and to conclude from an optically good fit that the kink-like behavior is absent. The 
kink-like feature might no longer be easily detectable by simple inspection of x systems 
where the entropy production approaches a nearly Gaussian form, but it still reveals itself 
in a maximum in 'x"{cf)- 

From our study, as well as from the observation of kink-like features in other systems 
fulfilling the fiuctuation theorem, we can therefore conclude that this behavior is generic for 
non-equilibrium systems. 



V. CONCLUSION 



Our study of interacting many-body systems has yielded two main results regarding the 
entropy production in the steady states of non-equilibrium systems. 

On the one hand, we showed that the kink observed in the entropy production large 
deviation function at zero entropy production is generic and directly follows from the fiuctu- 
ation theorem. The well known relation for the large deviation function ( ITSj) implies certain 
relations between the derivatives on both sides of a = 0, see Eq. f[T^ and flTB]) . It follows 
that the kink is characterized by the value of the mean entropy production rate in the steady 
state. 

On the other hand, our results indicate that the mean entropy production rate is a good 
measure for quantifying the non-equilibrium character of a system (or, equivalently, the 
degree of violation of detailed balance). Not only does it characterize the kink-like structure 
in the large deviation function, it also directly reveals whether a system gets closer to be 
an equilibrium system when changing some system parameters, as illustrated in Fig. [3^. 
In addition, the mean entropy production rate also reveals the presence of dynamic phase 
transitions, as seen in our study of the PASEP, see Fig. [3]d. 



Recently 



35H37| it was proposed that the non-equilibrium character of a system can be 



captured quantitatively by some global observables derived from the probability currents. 
Probability currents and entropy production are of course intimately related. They are 
unique to a non-equilibrium system and are absent in equilibrium systems. They also show 
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a remarkably similar dependence on system parameters. As an example, the reader should 



compare our Fig. [3^ with Fig. 2 in 26| where these quantities are computed for the one- 



dimensional reaction-diffusion systems discussed in this manuscript. To our knowledge, 
however, it has not yet been studied whether the observable derived from the stationary 
probability currents allows to probe the system for the presence of a phase transition. A 
more systematic comparison of the behavior of that quantity with that of the mean entropy 
production rate seems to be needed in order to further elucidate the commonalities and 
differences between these two quantities. 

As a final remark we note that in the past quite some effort has been devoted to the 
computation and understanding of the current large deviation function in driven systems. 
In our study we focused on the large deviation function of the entropy production and 
characterized its properties both for systems with (PASEP) and without (reaction-diffusion 
systems) particle currents. The exact relation between these two large deviation functions 
in driven systems remains to be better understood. 
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